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ESTIMATING APPARATUS AND METHOD OF INPUT AND OUTPUT 
ENABLING POWERS FOR SECONDARY CELL 



BACKGROUND OF THE INVENTION: 
5 Field of the invention 

[0001] The present invention relates to a 

technique to estimate powers which are enabled to be 
inputted to a secondary cell and which are enabled to 
be outputted from the same secondary cell. 

10 Description of the related art 

[0002] A Japanese Patent Application First 
Publication No. Heisei 9-171063 published on June 30, 
1997 exemplifies a previously proposed battery power 
calculating apparatus. In the previously proposed 

15 battery power calculating apparatus described in the 
above-identified Japanese Patent Application First 
Publication, an equation (V = R x I + Vo ) expressing 
an I-V straight line characteristic representing a 
discharge characteristic of the cell is calculated on 

20 the basis of a current I and a terminal voltage V 

supplied from a cell, an internal resistance R of the 
cell is calculated from its gradient, and an 
electromotive force Vo (which corresponds to a 
terminal voltage during a current interruption and 

25 also called an open voltage or open-circuit voltage) 
of the cell is calculated from an intercept. A 
minimum guarantee voltage value V m±n to guarantee a 
cell life on the basis of current I and cell 
temperature T is calculated and is substituted into 

30 the equation of I-V straight line to determine a 

maximum current value I ma x- The output enabling power 
value P is calculated from an equation of P = V ra m x 

^raax • 
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SUMMARY OF THE INVENTION: 

[0003] However, each of internal resistance R and 

open-circuit voltage V has a feature (characteristic) 
that each thereof R and V varies instantaneously (or 
5 continuously with respect to time) during charge-and- 
discharge operations in accordance with current I. In 
the above - described previously proposed power 
calculating apparatus disclosed in the above- 
identified Japanese Patent Application First 

10 Publication, current I and terminal voltage V are 
measured between two points during the charge 
operation in accordance with current I to calculate 
the I-V straight line. There is an assumption that 
internal resistance R and open-circuit voltage Vo 

15 determined from I-V straight line is not varied 
between two points. However, actually, since 
internal resistance R and open-circuit voltage Vo are 
instantaneously (or continuously) varied with respect 
to time, in the case of the calculation method 

20 disclosed in the above - described Japanese Patent 

Application First Publication, an estimation accuracy 
of output enabling power value P becomes lowered. 
[0005] It is, therefore, an object of the present 

invention to provide estimating apparatus and method 

25 for the secondary cell which are capable of 

estimating the input and output enabling powers for 
the secondary cell with a high accuracy and which are 
well (sufficiently) correspondent to an actual 
characteristic of the secondary cell. It is noted 

30 that the output enabling power is defined as a power 
which can be outputted from the secondary cell and 
the input enabling power is defined as the power 
which can be inputted into the secondary cell. 
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[0007] According to one aspect of the present 

invention, there is provided an estimating apparatus 
for a secondary cell, comprising: a current detecting 
section that detects a current (I) charged into and 
5 discharged from the secondary cell; a voltage 

detecting section that detects a terminal voltage (V) 
across the secondary cell; a parameter estimating 
section that integrally estimates all parameters (8) 
at one time in at least one of the following 

10 equations (1) and (2) with the measured current (I) 
and terminal voltage (V) inputted into an adaptive 
digital filter using a cell model described in a 
corresponding one of the following equations (1) and 
(2) whose parameters are estimated; an open-circuit 

15 voltage calculating section that calculates an open- 
circuit voltage (Vo) using the current (I), the 
terminal voltage (V), and the parameter estimated 
values (6); an input enabling power estimating 
section that estimates an input enabling power (P in ) 

20 of the secondary cell on the basis of the parameter 
estimated values (0) and open-circuit voltage (Vo); 
and an output enabling power estimating section that 
estimates an output enabling power (Pout) of the 
secondary cell on the basis of the parameter 

25 estimated values and the open-circuit voltage (Vo), 

the equation (1) being V = I +-^r7* Vo — U>. 

A(s) C(s) 

n n n 

wherein A(s) = Ya k • s k , B(s) = ^b k *s *. C(s) = ^c k * s ' 

Ic=Q k=Q *-0 

s denotes a Laplace transform operator, A(s), B(s), 
and C(s) denote each poly-nominal of s (n denotes 
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degrees) , a x * 0 , bi *0 , and c x ^ 0 and the equation 

(2) being V = +-2—*Vo ... (2), wherein 

A(s) A(s) 

n n 
A(s) = and B < s > = ' 

*-0 *=° 
[0008] According to another aspect of the present 

5 invention, there is provided an estimating method for 

a secondary cell, comprising: detecting a current (I) 

charged into and discharged from the secondary cell; 

detecting a terminal voltage (V) across the secondary 

cell; integrally estimating all parameters (6) at one 

10 time in at least one of the following equations (1) 
and (2) with the measured current (I) and terminal 
voltage (V) inputted into an adaptive digital filter 
using a cell model described in a corresponding one 
of the following equations (1) and (2) whose 

15 parameters are estimated; calculating an open-circuit 
voltage (Vo) using the current (I), the terminal 
voltage (V), and the parameter estimated values (8); 
estimating an input enabling power (P in ) of the 
secondary cell on the basis of the parameter 

20 estimated values (6) and open-circuit voltage (Vo); 

and estimating an output enabling power (P Q ut) of the 
secondary cell on the basis of the parameter 
estimated values and the open-circuit voltage (Vo), 

the equation (1) being V = ^7^ #/+ 7^T #Ko — t 1 *' 

A(s) C(s) 

n n n 

25 wherein A(s) = \a k 's k . B(s) = ^b k »s k . C(s) = ^c k »s . 

k=0 k=0 k=0 

s denotes a Laplace transform operator, A(s), B(s), 
and C(s) denote each poly-nominal of s (n denotes 
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degrees), ai * 0 , bi *0, and c x * 0 and the equation 

(2) being V = ^T #/+ ^7T #I ^ — < 2 )' wherein 

A(s) A(s) 

n n 
A(s) = Yak** and B(s) = ^b k •s . 

*= 0 *=0 

[0009] This summary of the invention does not 

5 necessarily describe all necessary features so that 
the invention may also be a sub - combination of these 
described features . 
BRIEF DESCRIPTION OF THE DRAWINGS: 

[0010] Fig.l is a functional block diagram of 

10 input and output enabling power estimating apparatus 
for a secondary cell according to the present 
invention applicable to each of first and second 
preferred embodiments . 

[0011] Fig. 2 is a specific circuit block diagram 

15 of a battery controller and a secondary cell load 

drive system to which the input and output enabling 
power estimating apparatus according to the present 
invention is applicable. 

[0012] Fig. 3 is a map representing a relationship 
20 between an open-circuit voltage and a charge rate 
( SOC ) . 

[0013] Fig. 4 is a model view representing an 
equivalent circuit model of the secondary cell in the 
input and output enabling power estimating apparatus 

25 of the first preferred embodiment. 

[0014] Fig. 5 is a model view representing an 
equivalent circuit model of the secondary cell in the 
input and output enabling power estimating apparatus 
of the second preferred embodiment. 

30 [0015] Fig. 6 is a processing flowchart 

representing a calculation process in the case of the 
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first preferred embodiment of the input and output 
enabling power estimating apparatus according to the 
present invention . 

[0016] Fig. 7 is a processing flowchart 
5 representing a calculation process in the case of the 
second preferred embodiment of the input and output 
enabling power estimating apparatus according to the 
present invention. 

[0017] Figs. 8A, 8B, 8C, 8D # 8E, 8F, 8G, 8H, 81 
10 show integrally a timing chart representing a result 
of a simulation based on the first embodiment of the 
input and output enabling power estimating apparatus 
according to the present invention. 

[0018] Figs. 9A, 9B, 9C, 9D, 9E, 9F, 9G, 9H, 91, 9J 
15 show integrally a timing chart representing a result 

of a simulation based on the second embodiment of the 

input and output enabling power estimating apparatus 

according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS: 
20 [0019] Reference will hereinafter be made to the 

drawings in order to facilitate a better 

understanding of the present invention. 

[0020] Fig. 1 shows a functional block diagram of 

input and output enabling power estimating apparatus 

25 according to the present invention for explaining a 

general concept of each of first and second preferred 
embodiments which will be described later. In Fig. 1, 
a reference numeral 1 denotes a parameter 6(k) 
estimating section that integrally estimates each 

30 parameter (the detailed description thereof will 

herein be omitted) in a cell model in which an open- 
circuit voltage Vo ( k ) is an offset term using 
measured voltage V and current I detected by current 
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I(k) detecting section 5 and terminal voltage V(k) 
detecting section 6. A reference numeral 2 denotes an 
open-circuit voltage calculating section Vo(k). The 
open-circuit voltage Vo(k) is calculated on the basis 
5 of the measured voltage V and current I and each 

estimated parameter. A reference numeral 3 denotes an 
input enabling power estimating section which 
estimates a power which can be inputted to the 
secondary cell on the basis of parameter 6(k) and 

10 open-circuit voltage Vo(k). A reference numeral 4 
denotes an output enabling power estimating section 
that estimates the power which can be outputted from 
the secondary cell on the basis of parameter 6(k) and 
open-circuit voltage Vo(k). A reference numeral 5 

15 denotes a current I(k) detecting section that detects 
the current charged to or discharged from the 
secondary cell. A reference numeral 6 denotes a 
terminal voltage V(k) detecting section that detects 
a terminal voltage of the cell. 

20 [0021] Fig. 2 shows a block diagram representing 

a specific structure of a battery controller and a 
secondary cell load driving system to which the input 
and output enabling power estimating apparatus 
according to the present invention is applicable. In 

25 this system, the input and output enabling power 

estimating apparatus is mounted in a system in which 
a load such as a motor is driven and a regenerative 
power of the motor is used to charge the secondary 
cell. In Fig. 2, a reference numeral 10 denotes a 

30 secondary cell (or merely, a cell), a reference 

numeral 20 denotes a load of the motor or so on, and 
a reference numeral 30 denotes a battery controller 
(an electronic control unit) which functions to 
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estimate the input and output enabling powers of cell 
10. Battery controller 30 includes a microcomputer 
including a CPU ( Central Processing Unit) that 
calculates a program, a ROM (Read Only Memory) that 
5 stores the program, a RAM (Random Access Memory) 
storing a result of calculations, and electronic 
circuits. A reference numeral 40 denotes a current . 
meter that measures (detects) a current which is 
charged into and discharged from secondary cell 10. A 

10 reference numeral 50 denotes a voltage meter to 

detect a terminal voltage across secondary cell 10. 
These meters are connected to battery controller 30. 
The above-described battery controller 30 corresponds 
to parameter 0(k) estimating section 1 of Fig. 1, 

15 open-circuit voltage Vo(k) calculating section 2, 

input enabling power estimating section 3, and output 
enabling power estimating section 4. In addition, 
current meter 40 corresponds to current I ( k ) 
detecting section 5 and voltage meter 50 corresponds 

20 to a terminal voltage V(k) detecting section 6, 

respectively. It is noted that a reference numeral 60 
shown in Fig. 2 denotes a temperature sensor to 
detect a cell temperature and a reference numeral 70 
shown in Fig. 2 denotes a relay circuit (or simply a 

25 relay) . 

[0022] (First Embodiment) 

Next, a, so-called, cell model used in the 
first embodiment will be described below. Fig. 4 
shows an equivalent circuit model of the secondary 

30 cell in the first embodiment. This equivalent circuit 
model corresponds to a case where denominators of 
right side first term and right side second term are 
the same as shown in equation (2). This equivalent 
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circuit model is a reduction model (first order or 
fist degree) in which positive pole and negative pole 
are not especially separated) but enabled to 
represent a relatively accurate charge - and- dis charge 
5 characteristic of the actual cell. In Fig. 4, a 
model input is current I [A] (Amperes) (positive 
value is charge and negative value is discharge) and 
model output is a terminal voltage V[V]. Vo[V] 
denotes an open-circuit voltage (or called, an 
10 electromotive force or open voltage). K denotes an 

internal resistance, Ti and T 2 denote time constants. 
The cell model can be expressed in the following 
equation (3). It is noted that s denotes a Laplace 
transform operator. 



15 



v m + —k-.y 0 ... (2), wherein 

A(s) A(s) 

A(s) "2 fl * #5 * ' B(S) = ti 1 **'** 

It is noted that A(s) and B(s) denote polynominals of 
s, n denotes a degree (order number), and ai *0 
and bi *0 . 

20 V = — + - 7 # ^o — ( 3 )- 

T^s+l ^'5 + 1 

It is noted that equation (3) is a variation of 

equation (2) in which Ti*s + 1 is substituted for 

and A(s) (A(s) = T x * +1)) and K •(T 2 *s + 1) is 

substituted for B(s) (B(s)= K •(T 2 *s + 1)). In the 

25 case of such a lithium ion battery that a 

convergence of the open-circuit voltage is relatively 

fast, the denominators of right side first term and 

right side second term can be represented by the same 

time constant Ti as appreciated from equation (3). 
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[0023] Hereinafter, a procedure of a derivation 

from the cell model in equation (3) to an adaptive 
digital filter will first be explained. The open- 
circuit voltage Vo can be written in the following 
5 equation (4) assuming that a value of current I 

multiplied by a variable efficiency h is considered 
to be an integration value from a certain initial 
state, 

Vo = -•/ ...(4) . 

s 

10 If equation (4) is substituted into equation (3), the 
following equation (5) is given. If equation (5) is 
rearranged, the following equation (6) is given. If 
a stable low pass filter Gi p (s) is multiplied to both 
sides of equation (6) and rearranged, then, the 

15 following equation (7) is given. 



K*(T 2 *s+l) T 1 h _ 

V = — -•/ + • — •/ ... (5). 

^•5+1 s 

K*T 2 »s 2 +K*s+h h T 

V = 1 • — •/ ... (6). 

T x • s + 1 s 

Gi p (s)*(Ti'S 2 + s)*V = G lp ( S ) • (K*T 2 *s 2 + K-s + h)*I 
20 ... ( 7 ) . 

A value of an actually measurable current I and 
terminal voltage V for which a low pas filter and a 
band pass filter are processed is defined as in the 
following equation (8). A time constant p of equation 
25 (8) is a constant determining a response 
characteristic of 6i p (s). 
[0024] 

V 2 -s**G ¥ i*yV V 2 =s-G lp (s)-V Vt-G+WV G ip = {p .l +1 y 

I>=s 2 *G lp (s)'I I 2 =s-G lp (s)-I I, =G Ip (s)-I 
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... (8), 

[0025] If, using equation (8), equation (7) is 

rewritten, then, the following equation (9) is given. 
Furthermore, if equation (9) is deformed, an equation 
5 (10) is given. 

T x -K 3 +K 2 =K'T 2 •I 3 +K*I 2 +h*I 1 ... ( 9 )• . 

V 2 =-r i *K 3 +^T 2 •I,+K*I 2 +d*I 1 = [V 3 I 3 I 2 /J- 

...(10). 

Equation (10) indicates a product - and- sum equation 
10 between measurable value and unknown parameters. 

Hence, equation (10) is coincident with a standard 
form (equation (11)) of a general adaptive digital 
filter. It is noted that, in equation (11), y = V 2 , 
u T = h h 'iL 9 T = t-Ti K • T 2 K h]. 

15 y = (0 T • 0 ...(H) . 

[0026] Hence, if current I and terminal voltage V 

to both of which a filter is processed are used for 
the adaptive digital filter calculation, unknown 
parameter vector 0 can be estimated. In this 

20 embodiment, a, so-called, both eyes trace gain method 
which has improved a logical defect of the adaptive 
filter (namely, once an estimated value is converged, 
an accurate estimation, thereafter, cannot be made 
any more even if the parameter is varied) is used. 

25 [0027] Upon an assumption of equation (11), a 

parameter estimation algorithm to estimate an unknown 
parameter vector 0 is described in an equation (12). 
It is noted that parameter estimated values at a time 
point of k is assumed to be 0(k). 



K*T 2 

K 
h 
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y(k) = 



A 3 



1 + A 3 •o) T (k)^P{k-\)»co{k) 
6(k) = 6(k -1) - y (k) • P(* - 1) • co(k) •[co T (k)* d(k -1) - y(k)] 



P(k) = 



P(k-1)- 



X 3 'P(k-l)»co(k)» co (k) * P(k - 1) 



1 + A 3 • o) r (k) • P(k -1) • <u(A:) 



\(k) = 



frflcejp'CA:)} . A ^ trace\p\k)} 
YU YU 



trace 



YU 



trace 



H*)} . trace\p\k)} 



YL 



Yl 



Yl 



... (12). 

In equation (12) , Xi , X 3 , Yu , Yl denotes initial set 
values, 0 < Xi < 1, 0 < X 3 < °°. In addition, P(0) is 
5 a sufficiently large initial value and tracejp} means 
a trace of a matrix P. In this way, the derivation of 
the adaptive digital filter from the cell model has 
been explained. 

[0028] Fig. 6 shows an operational flowchart 

10 executing a microcomputer of battery controller 30. 
The routine shown in Fig. 6 is repeated for each 
constant period of time T 0 . For example, I(k) is the 
present value and I(k-l) is the value one time before 
the present time of k. 
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At a step S10A, battery controller 30 measures 
current I(k) and terminal voltage V(k). At a step 
S20A, battery controller 30 determines whether an 
interruption relay of secondary cell 10 is in an 
5 engaged state (closed) or in an interrupted state 

(open). It is noted that battery controller 30 also 
controls the interruption relay. If the relay is 
interrupted (current 1=0), the routine goes to a 
step S30A. If the relay is engaged, the routine goes 

10 to a step S40A. At step S30A, battery controller 30 
determines that terminal voltage V(k) is stored as 
terminal voltage initial value V_ in t . At step S40A, 
battery controller 30 a difference value of the 
terminal voltage AV(k): A V ( k ) = V(k) - V_ ini 

15 Since this the initial value of the estimation 

parameter within the adaptive digital filter is set 
to about zero, the inputs are all zeroed to prevent 
each estimated parameter during a start of an 
estimation calculation from being diverged. Whenever 

20 the relay is interrupted, the routine goes to step 
S30A. Hence, since 1=0 and AV(k) = 0, the 
estimated parameters still remains in the initial 
state. At a step S50A, battery controller 30 
performs a low pass filter (LPF) and a band pass 

25 filter (BPF) processing based on equation (13) for 
current I(k) and terminal voltage difference value 
AV(k) to calculate Ii through I 3 and Vi through V 3 . 
At this time, in order to improve an estimation 
accuracy of the parameter estimation algorithm of 

30 equation (12), a response characteristic of low pass 
filter Gip(s) is set to be slow so as to reduce 
observed noises. 
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[ 0029 ] 



It is noted that this response 



characteristic is made faster than the response 
characteristic of the cell. Time constant p in 
equation (13) is a constant to determine the response 
5 characteristic of Gi p (s) during the equation (13). 



10 S50A into equation (12) to calculate parameter 

estimated value 8(k). It is noted that y = V 2 , co T 
= [V 3 I 3 I 2 111. and 0 T = [-Ti K*T 2 K h]. 
At a step S70A, back-up controller 30 substitutes T 1 , 
K • T 2 , and K from among parameter estimated value 

15 6(k) calculated at step S60A into the following 

equation (14) from among parameter estimated values 
6(k) calculated at step S60A into equation (14), Ii, 
and I 2 , and Vi , and V 2 calculated at equation (13). 
V 0 =(T X •5+l)*K-iC-(r2*5 + l)*/ 
AV 0 =G lp (s) 

= g 1p (s)*{(ti*s+i)*v-k*(t 2 *s+i)*i} 

20 ... (14). 

Equation (14) is a deformation for the cell model 
(equation (3)) and low pass filter G ip (s) is 
multiplied to both sides. Then, voltage component of 
AVo is replaced with open-circuit voltage Vo (Vo is 
25 substituted for AVo ) . Since the variation of open- 




V 2 =s*G lp (s)'V, V 1 =G lp (s)*V 
I 2 =s»G lp (s)*I, h=G lp (s)*I 



h -s 2 -G lp (s)^I, 



...(13) . 

[0030] At a step S60A, controller 30 substitutes 

Ix through I 3 and Vi through V 3 calculated at step 
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circuit voltage Vo is moderate, Vo can be replaced as 
follows: AVo = G ip (s) # Vo. 

Since variation rate AVo(k) of the open-circuit 
voltage estimated value from a time at which the 
5 start of the estimated calculation is carried out, 
the initial value at the later stage of a step S80A. 
[0031] At step S80A, the open-circuit voltage 

initial value, namely, a terminal voltage initial 
value V_ lni is added to AVc(k) calculated at step 

10 S70A to calculate open-circuit voltage estimated 
value Vo(k) using the following equation (15). 

V 0 (k) = AV 0 (k) + V_ lri i ... (15). 

At a step S90A # battery controller 30 calculates a 
charge rate SOC(k) from Vo(k) calculated at step S80A 

15 using a correlative map of the open-circuit voltage 

shown in Fig. 3 and charge rate. It is noted that VL 
shown in Fig. 3 is an open-circuit voltage 
corresponding to SOC = 100 % and VH shown in Fig. 3 
is an open-circuit voltage corresponding to SOC = 

20 100 %. 

[0032] At a step S100A, battery controller 30 

calculates an input enabling power estimated value 
P ln and an output enabling power estimated value P Q ut • 
Hereinafter, the detailed description of the 

25 calculation method of the input enabling power 
estimated value will be described below. 

In the cell model (equation (3)), in a case 
where a transient characteristic is ignored, equation 
(16) is resulted. This means that this means a 

30 quantitative cell model. 

V = K • I + Vo ... ( 16 ) . 
[0033] Suppose that the terminal voltage of the 

cell immediately before a predefined excessive (or 
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over) charge is resulted is a maximum enabling 
voltage V ma x and the terminal voltage of the cell 
immediately before the predefined excessive (or over) 
discharge is resulted in a minimum enabling voltage 
5 V min . Then, in order to calculate the input enabling 
power estimated value Pin, it is necessary to 
require the current value by which the terminal 
voltage has reached to maximum enabling voltage V max - 
Hence, using equation (16) in which the transient 
10 characteristic is ignored and equation (16) is used 
to calculate maximum input current Im_max using 
equation (16). 

[0034] In equation (16), maximum enabling voltage 

V max is substituted into V, estimated value K from 

15 among the parameter estimated values 6(k) calculated 
at step S60A is substituted into K and circuit 
voltage estimated value Vo(k) calculated at step S80A 
is substituted into Vo , respectively, to calculate a 
maximum input current I ln_max • 

20 [0035] In the same way as the case of output 

enabling power estimated value Pout , minimum enabling 
voltage V min is substituted into V in equation (16), 
estimated value K from among the parameter estimated 
value 6(k) calculated at step S60A is substituted 

25 into K, and circuit voltage estimated value Vo ( k ) 
calculated at step S80A is substituted into Vo , 
respectively, to calculate maximum input current 
Im_max. Then, input enabling power estimated value 
Pin and output enabling power estimated value P ou t are 

30 calculated from equation (17), 
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Pin = I *V 

* * * m _ max max 

V_„-Vo 



=1/ 



<?u/ max 



9 m 



min 



min 



...( 17 ) . 

Maximum enabling voltage V max is a terminal voltage in 
a case where the cell is charged to a voltage 
5 immediately before the cell is the excessive charge. 
Minimum enabling voltage V min is a terminal voltage in 
a case where the cell is discharged to a value 
immediately before the cell is the excessive charge. 
These maximum enabling voltage V ma x and minimum 

10 enabling voltage V m ± n are variables determined by the 
kind of cells and the cell temperature. For example, 
a relationship between the cell temperature and V max 
determined according to, for example, the experiments 
and a relationship between the cell temperature and 

15 V min can be stored as maps and a map reference can be 
used to calculate V max and V min . At a step S110A, 
numerical values required for the subsequent 
calculation are stored and the present calculation is 
ended. An operation of the first embodiment has been 

20 described above. 

[0036] Hereinafter, an action and advantages of 

the estimating apparatus for the secondary cell in 
the first embodiment will be described below. 
[0037] In the first embodiment, since the 

25 relationship between current I of the secondary cell, 
terminal voltage V, and the open-circuit voltage Vo 
approximates the transfer function such as in 
equation (2). Specifically, equation (3), it becomes 
possible to apply the adaptive digital filter (well 
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known estimating algorithm) such as the method of the 
least square. Consequently, it becomes possible to 
integrally estimate the parameters in equations 
(coefficients of poly-nominals (A(s) and B(s)). When 
5 the estimated parameters are substituted into 

equation (2), the estimated value of open-circuit 
voltage Vo can easily be calculated, 

[0038] These unknown parameters are affected by 

the charge rate (SOC), the cell temperature, and a 

10 degree of deterioration. Although these parameters 
are known to be instantaneously varied with respect 
to time, the adaptive digital filter can sequentially 
be estimated with a high accuracy. Since input 
enabling power P in and output enabling power P out are 

15 estimated using the estimated coefficient parameters 
and the open-circuit voltages Vo, the input and 
output enabling power P in and P out can be estimated, 
even if the input and output enabling powers are 
varied during the charge or discharge operation, its 

20 variation can accurately follow to estimate the input 
and output enabling powers . 

[0039] As compared with the second preferred 

embodiment as will be described later, since an 
easier cell model ( equat ions ( 2 ) and (3)) is used, a 

25 formalization (or an equalization) of the adaptive 

digital filter becomes easy and the numbers of times 
the calculations are carried out can be reduced. 
[0040] Figs. 8A through 81 show integrally results 

of simulations of the input and output enabling power 

30 estimations based on the first embodiment. 

[0041] In Figs. 8A through 81, with a time of 400 

seconds as a boundary, the cell parameters are 
changed in a stepwise manner from a high temperature 
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corresponding value to a low temperature 
corresponding value. It is noted that, in this 
example of Figs. 8A through 81, such as a lithium ion 
battery, the cell having a fast convergence of the 
5 open-circuit voltage is presumed. As appreciated 

from Figs. 8A through 81, time constants T if T 2 , and 
internal resistance K are considerably coincident 
with real values even if the cell parameter given 
when the simulations are carried out are varied in a 

10 stepwise manner. Hence, the open-circuit voltage 
estimated values are similarly coincident with the 
true values. In the first embodiment, using the 
estimated coefficient parameter, the open voltage Vo , 
and maximum enabling voltage V m in, input enabling 

15 power P in is estimated. Hence, even if cell 
parameters and open-circuit voltage Vo are 
instantaneously varied with respect to time during 
the charge or discharge operation, the output 
enabling power estimated value is accurately 

20 coincident with the true value (real value). 
[0042] (Second Embodiment) 

Next, an operation of the second preferred 

s 

embodiment will be described. First, the cell model 
used in the second embodiment will be explained below. 
25 Fig. 5 shows an equivalent circuit model of the 
secondary cell in the second embodiment . 
[0043] Before explaining the equivalent circuit 

model shown in Fig. 1, equation (1) is described 
herein . 

30 v^BQ.^^—.Vo ... (1), 

A(s) C(s) 
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wherein A(s), B(s), and C(s) denote a poly-nominal of 
s (n denotes an order number) and sl% * 0 , bi*0 , and ci 
* 0 . 

The equivalent circuit model corresponds to a case 
5 where the denominators of the first term and the 
second term are different as described in equation 
(1). This equivalent circuit model is a reduction 
model (first degree or first order) in which a 
positive pole and a negative pole are not specially 

10 separated from each other but can relatively 
accurately indicate the charge-and-discharge 
characteristics of the actual cell. In Fig. 5, a 
model input is current I [A] (positive value 
indicates the charge and negative value indicates the 

15 discharge) and a model output indicates terminal 
voltage V [V] and Vo[V] indicates open-circuit 
voltage (referred also to as an electromotive force 
or open (circuit) voltage). A symbol K denotes the 
internal resistance. Ti through T 3 denotes time 

20 constants. This cell model can be expressed as in the 
following equation (18). It is noted that s denotes 
a Laplace transform operator. As a lead-acid battery 
cell ( or lead storage battery) , the convergence of 
the open-circuit voltage is very slow cell, the 

25 relationship between Ti « T 3 is present. 

K .*-(T a -5 + l). / + _l 

^•5+1 r 3 *5 + i 

Equation (18) is a replacement of equation (1) with 
A(s) = Ti*s + 1, B(s) = K»(T 2 # s + 1). First, the 
derivation of the cell model shown in equation (18) 
30 to the adaptive digital filter will be explained 
below . 
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The open-circuit voltage Vo can be written with the 
value of current I multiplied by a variable 
efficiency h considered as an integration value from 
a certain initial state, 

5 V 0 =-•/ (19)- 

s 

If equation (19) is substituted into equation (18), 
the following equation (20) is given. If arranged, 
the following equation (21) is given. 

V = — -•/ + • — •/ ... (20). 

^•5+1 T 3 's+l s 

io 5*(rj •5 + i)(73 •s+i)*v = K'(T 2 •5+i)(r 3 •5+i)*5-/+/i*(r 1 

T x T 3 •s 3 +(T X + T 3 )*s 2 + = 

K>r 2 T 3 + A>(T 2 + r 3 )*s 2 + (K + h*T 1 )* s + ^K + h'T^* s+h\'I 
( a»s 3 +b*s 2 +s ) *V = (c»s 3 + d»s 2 +e # s + f) # I ...(21). 
It is noted that the parameters shown in equation 

15 (21) is rewritten as follows: 

a = T 1* T 3* b = Ti + T 3 , c = K*T 2 *T 3 

d - K*(T 2 + T 3 ), e = K + h*T lr f = h ...(22). 

If a stable low pass filter Gl(s) is introduced into 

both sides of equation (21) and arranged, the 

20 following equation (23) is given. 

1 (a*s 3 +b*s 2 +s)'V = (c*5 3 +d • s 2 +e's + f)*I ...(23) 



G^sy G(s) 
A value of each of actually measurable current I and 
terminal voltage V for which low pass filter is 
processed and band pass filter is processed is 
25 defined as shown in an equation (24). In equation 
(24), pi denotes a time constant determining the 
response characteristic of G x (s). 
[ 0044 ] 
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G 1 (s) 



•I 



1 G^s) 



'I 



G 1 (s) 



G x (5) 



G!( S ) 



G^s) 



Gi(5) 
... (24). 

If equation (23) is rewritten using variables shown 
in equation (24), an equation (26) is resulted. If 
deformed, the following equation (26) is given. 
a«K 3 +b'V 2 +Vi =c*I 3 +d • I 2 +e»I x + f • I 0 ; and 

Vi= -a»V 3 -b»V 2 + c»I 3 + d»I 2 + e # Ii + f*Io ...(25). 



- a 

-b 

c 

d 

e 

f 



...(26) 



Since equation (26) indicates a product - and- sum 
equation between measurable values and unknown 
parameters, equation (26) is coincident with a 
standard form (equation (27)) of a general adaptive 
digital filter. It is noted that (o T means a 
transposed vector in which a row and a column of a 
vector (o are replaced with each other. 

y = co T *0 ... {21) . 

It is noted that y = V ir 
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- a 

-b 
c 
d 

e 

f . 

Hence, using the adaptive digital filter calculation 
is used for the filter processed signals by which 
current I and terminal voltage V are filter processed 
so that the unknown parameter vector 9 can be 
estimated. In this embodiment, the simple, so-called, 
both eyes trace gain method is used which improves 
the logical defect (the accurate estimation cannot 
again be made once the estimated value is converged) 
of the adaptive filter by means of the least square 
method is used. On the premise of equation (27), a 
parameter estimation algorithm to estimate unknown 
parameter vector 0 is given as in the following 
equation (28). It is noted that parameter estimated 
value at a time point of k is assumed to be 9(k). 



co T = [v 3 V 2 I 3 I 2 h A)L e = 
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1 + A 3 • co T (k) • P(k -1) • cu(&) 



- 6(k - 1) - y (k) • P(A: - 1) • co(k) • \co (k) • 6{k - 1) - y (&) 



4 



l + A 3 »a) T (k)»P(k-l)'co(k) 



trace\ 



re/ 



trace{P'(k)} , trace 



YU 



trace 



M. *rflce{p'(*)} 



Yl 



YL 



Yl 



... ( 28 ) . 

It is noted that X lf X 3 , Yu , and Yl are initial set 
values and 0 < X ± < 1 and 0 < X 2 (k) < <* . P(0) has a 
5 sufficiently large initial value and has a 

sufficiently small initial value which is not zero. 
trace{P} means a matrix P trace. In this way, the 
derivation of the cell model to the adaptive digital 
filter has been described. 

10 Fig. 7 shows an operational flowchart carrying out a 
microcomputer of battery controller 30. The routine 
shown in Fig. 7 is executed whenever a constant 
period of time To has passed. For example, I(k) means 
the present value and I(k-l) denotes one previous 

15 value as described in the first embodiemnt . In Fig. 7, 



25 



10 



the contents of steps S10B through S40B are the same 
as those of steps S10A through S40A described in Fig 
6. Hence, the explanation thereof will herein be 
omitted. At a step S50B, the low pass filtering and 
the band pass filtering are carried for current I(k) 
and terminal voltage difference value AV ( k ) on the 
basis of the following equation (29) to calculate 
Io(k) through I 3 (k) and V L (k) through V 3 (k). 
[0045] 



(his) 

s 



<h<?) 



It. = 



Gj(5) 



V, = 



G 1 (s) 



V 0 = 



G 1 (s) 



'V 



K, =• 



... ( 29 ) . 

It is noted that, in this case, in order to improve 
an estimation accuracy of parameter estimation 

15 algorithm of equation (28), the response 

characteristic of low pass filter Gi(s) is set to be 
slow so as to reduce the observed noises. However, if 
the response characteristic of low pass filter d ( s ) 
is faster than the response characteristic (an 

20 approximate value of time constant Ti is already 

known) of the cell model, each parameter of the cell 
model cannot accurately be estimated. Pi in equation 
(29) is a constant determining response 
characteristic of Gi(s). 
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25 



[0046] At a step S60B, controller 30 substitutes 

I 0 (k) through I 3 (k) and Vi(k) through V 3 (k) into 
equation (28). The calculation in accordance with 
equation (28) which is the parameter estimation 
algorithm is, then, calculated to determine parameter 
estimated value 6(k). It is noted that y(k), co T ( k ) , 
and 0(k) are given in the following equation (30). 
[0047] 

co T (k) = [v 3 (k) V 2 (k) I 3 (k) I 2 (k) I x (k) I Q (k)] 



0{k) = 



-a(k)] 
-b(k) 
c(k) 
d(k) 
e(k) 

/(*) 

... ( 30 ) . 

At a step S70B, the filtering processes of low pass 
filter and band pass filter are carried out on the 
basis of current I(k) and terminal voltage difference 
value AV(k) on the basis of an equation (34) to 
calculate I 4 (k) through I 6 (k) and V 4 (k) through V 6 (k). 
a through e from among parameter estimated values 
9(k) calculated at step S60B are substituted into 
equation (33) which is a deformation from equation 
(18) to calculate AVo which is used in place of the 
open-circuit voltage Vo . Since the variation in 
open-circuit voltage Vo is moderate, AVo can be 
substituted. It is noted that the derivation at step 
S70B is the variation quantity AVo ( k ) of the open- 
circuit voltage Vo ( k ) from a time at which the 
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estimation calculation is started. Therefore, the 
initial value is added at a later step S90B. It is 
noted that, at the derivation of equation (33), K in 
the equation (32) and e of the equation (33) are 
5 strictly different from each other. Physically, K » 
h • Ti, e is approximated to K (e is about equal to K, 
e s K). In addition, since the approximate value of 
T r of the cell parameter is known as several seconds, 
ti in equation (34) is set to a value near to the 
10 approximate value of T L . Thus, since a term of (Ti • s 
+ 1) which is rested on a numerator in equation (33) 
can be cancelled, the estimation accuracy of open- 
circuit voltage Vo can be improved. 

15 • V Q = V ^ — - • I 

T 3 's + 1 u Ti-s + l 

(Ti •5+l)»K 0 =(r x '5+1)^3 *s + l)*V -K'(T 2 's+l)(T 3 
(7i • S + 1) • V 0 = { (Ti • T 3 • 5 2 + (Tx + T 3 ) • 5 + 1 } • V 

-lf»T 2 'T 3 's 2 + K»(T 2 +T 3 )*s + k\'I ... ( 31 ) . 



G 2 (s) ° G 2 ( S y G 2 (s) K 
...(32) . 

AV 0 = (ri ' Hl) >V Q = a'V 6 +b*V 5 +V 4 -c»I 6 -d'I 5 -e«/ 4 ...(33) 
G 2 (s) 

In equation (33), a = T x • T 3 , b = T x + T 3 , c = K 
•T 2 »T 3 , d = K »(T 2 + T 3 ), c = K + h # Ti = K. 
25 [0048] 
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4 C 2 (j) 4 G 2 (s) 



<?2(*) G 2 (5) G 2 (s) P2*s+1 n-j+l 

G 2 C0 G 2 (j) 
...(34) . 

When the calculated AVo(k) is substituted into an 
equation (35), estimated value AVo'(k) only at a 
5 right side second term of cell model (refer to 
equation (18)) is calculated, Vo(k) denotes an 
estimate value of the open-circuit voltage itself 
approximated by equation (18) and Vo' denotes an 
apparent estimated value of the open-circuit voltage 
10 appearing on the terminal voltage. It is, however, 

noted that, in the derivation of equation (35), T 3 at 
the right side is strictly different from right side 
b. Physically, since T 3 » ?i. b = T 3 + T x ^ T 3 . 

AJ/'o = • AV Q s *AV 0 ...( 35) . 

T 3 *s+1 ° fs+1 

15 Equation (35) corresponds to Vo/C(s). That is to say, 
Vo = AVo and C(s) = T 3 • s + 1 ^ b • s + 1 . 
[0049] At a step S80B of Fig, 7, the open-circuit 

voltage initial value, namely, terminal voltage 
initial value V_ in i is added to each of Vo ( k ) and 

20 Vo'(k) calculated at step S70B. That is to say, open- 
circuit voltage estimated value Vo ( k ) is calculated 
using the equation (36) and apparent open-circuit 
voltage estimated value Vo'(k) is calculated using 
the following equation (37). It is noted that 

25 estimated value Vo 1 is not the estimated value Vo 1 of 
open-circuit voltage Vo itself but the apparent open- 
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circuit voltage estimated value appearing on the 
terminal voltage. 

[0050] V 0 (k) = AVi(k) + V_ ini ...(36). 

Vo'(k) = AVo'(k) + V_ ini ... (37). 

5 At a step S90B, battery controller 30 calculates 

charge rate SOC(k) from Vo(k) calculated at step S80B 
using the correlation map between the open-circuit 
voltage and charge rate shown in Fig. 3. It is noted 
that V L shown in Fig. 3 is the open-circuit voltage 

10 corresponding to SOC = 0% and V H is the open-circuit 

voltage corresponding to SOC = 100%. At a step S100B, 
battery controller 30 determines whether estimated 
value Vo(k) is equal to or larger than open-circuit 
voltage estimated value Vo'(k). This determination at 

15 step S100B functions to search either one of which is 
nearer to either maximum enabling voltage V max or 
minimum enabling voltage V m m . It is noted that 
maximum enabling voltage V max or minimum enabling 
voltage V min is a variable determined by the kinds of 

20 cells and temperature in the cell. The calculation 
method is the well known method so that they can be 
determined using the well known technique as 
described in the first embodiment. At step S100B, if 
Vo'(k) ;> Vo(k), the routine goes to a step S110B. If 

25 Vo'(k) < Vo(k), the routine goes to a step S120B. At 
a step S110B, battery controller 30 calculates input 
enabling power estimated value P in and output 
enabling power estimated value Pout • In the cell model 
(equation (18), the cell model is expressed in 

30 equation (38) in a case where a transient 
characteristic is ignored and this means a 
quantitative cell model. To calculate input enabling 
power estimated value P± n , the current value to reach 
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to maximum enabling voltage V max is needed. Hence, 
maximum input current I in _max is calculated using 
equation (38) in which the transient characteristic 
is ignored. That is to say, at step S110B, since 
5 Vo'(k) a: Vo(k), Vo'(k) is nearer to maximum enabling 
power V m ax and Vo(k) is nearer to minimum enabling 
voltage V min • Hence, in order to calculate input 
enabling power estimated value P in , maximum enabling 
power voltage V ma x is substituted for V of equation 

10 (38), estimated value e from among parameter 

estimated values 9(k) calculated at step S60B is 
substituted for K of equation (38), and Vo'(k) 
calculated at step S80B is substituted for Vo of 
equation (38) so that maximum input current I± n _max is 

15 calculated from equation (38) obtained from equation 

(39) . 

[0050] V = K«I+V 0 ...(38). 

V max = e • Iin.max + V 0 ' ...(39). 

[0051] On the other hand, for output enabling 

20 power estimated value P ou t , minimum enabling voltage 

V min is substituted into V, one of parameter estimated 
values, viz., e from among parameter estimated 
values 0(k) calculated at step S60B is substituted 
into K, open-circuit voltage estimated value Vo(k) 
25 calculated at step S80B is substituted for Vo of 

equation (38). The obtained equation is an equation 

(40) to calculate maximum output current I 0 ut_max- 

V m in = e * Iout_max + V Q ... (40). 

[0052] Next, using maximum input current Im_max. 

30 maximum output current I 0 ut_max derived as described 
above, equations (41A and 41B) calculate an input 
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enabling power estimated value (P in ) and an output 
enabling power estimated value ( P out ) . 

It is noted that, at the derivation of maximum input 
current I in _ max and maximum output current I 0 ut_max, K in 
equation (38), e in equations (39) and (40) are 
strictly different from one another. However, since, 
physically, K » h • Ti, e = K + h*T x ^K. 

P- = /• •V 

J- ift * tn ttiqy r rr 



1 in max r max 



^max . t/ 



max 



max 



max 



...(41A). 



Pout -I ^ 



out max 



min 



V -V • 
o Y min # 



min 



.(41B). 



10 [0053] At a step 120B, battery controller 30 

calculates input enabling power estimated value P in 
and output enabling power estimated value P ou t • Since 
step S120B is the case where Vo'(k) < Vo(k), Vo(k) is 
nearer to maximum enabling voltage V max and Vo'(k) is 

15 nearer to minimum enabling voltage V min . Hence, in 
order to calculate input enabling power estimated 
value P in , maximum enabling voltage V max # estimated 
value e from among parameter estimated value 0(k) 
calculated at step S60B using an equation (42) 

20 obtained by substituting Vo(k) calculated at step 
S80B into equation (38). Thus, an equation (43) is 
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given. Maximum output enabling current I 0 ut_max is 
calculated using equation (43). 

[0054] Vmzx = e'I in max +K D ... (42). 



^min ~ e * I* 



out max + 



(43) 



Next, using maximum input current Im.max and maximum 
output current I 0 ut_max, input enabling power estimated 
value P in and output enabling power estimated value 
P out are calculated from equations (44A) and (44B) as 
will be described below. 



10 [0055] 



I. •V 

in _ max max 

V -V 

max 0 m rs 



( 44A) 



Pout _ l lout max I # ^min 



V -V • 

Y o v min # t/ 



min 



— — --v, 



b 9 s+l 



mm 



( 44B) 



At a step S130B, battery controller 30 stores 
numerical values needed for the next calculation and 

15 the present calculation is ended. The second 

preferred embodiment of the estimating apparatus 
according to the present invention has been described 
[0056] Next, the action and advantages of the 
second embodiment of the estimating apparatus will be 

20 described below. In the second embodiment, the 

relationship from among current I of the secondary 
cell, terminal voltage V, and open-circuit voltage Vo 
is constituted to be approximated by means of the 



-33- 



transfer function such as equation (1) (specifically, 
equation (18)), it is possible to apply to the 
adaptive digital filter of the method of least 
squares. Consequently, it becomes possible to 
5 integrally estimate the parameters at one time 

(coefficients of poly - nominals of A(s), B(s), and 
C(s)). Since the estimated parameters are 
substituted into equation (1), the estimated value of 
open-circuit voltage Vo can easily be calculated. 

10 These unknown parameters are affected by a charge 
rate (SOC, viz.. State Of Charge), an ambient 
temperature of secondary cell, and a degree of 
deterioration is varied instantaneously 
(continuously) with respect to time. However, the 

15 sequential estimation can be made with a high 

accuracy by means of the adaptive digital filter. 
Since input enabling power P in and output enabling 
power Pout are estimated using the estimated 
coefficients (parameters) and the open-circuit 

20 voltage Vo . Hence, even if, together with the 

variation in cell parameters during the charge-and- 
discharge operation, the input and output enabling 
powers P in and P ou t are varied, the adaptive digital 
filter can accurately follow its variation so that 

25 the input and output enabling powers can accurately 
be estimated. 

[0057] Fig. 9A through 9J integrally show a 

simulation result of the input output enabling power 
estimations based on the second embodiment. In Fig. 
30 9A through 9J, with a time of 500 seconds as a 

boundary, the cell parameters are varied from a low 
temperature corresponding value to a high temperature 
corresponding value in a stepwise manner. In the 
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case of the simulation, as far as a time constant of 
a first-order lag described in the cell model 
(equation (18) is concerned, T x « T 3 is set. This is 
because, the cell having the very slow convergence 
5 characteristic of open-circuit voltage Vo like the 
lead-acid battery is assumed and set. 

[0058] As appreciated from Figs. 9A through 9J, 

parameter estimated values a through e that the 
adaptive digital filter outputs are coincident with 

10 their real values even if the cell parameter given 
when the simulation is carried out is varied in the 
stepwise manner (substantially at a right angle). 
Hence, the open-circuit voltage estimated value are 
coincident with the real values. In the second 

15 embodiment, input enabling power P in is estimated 
using the estimated coefficient parameter, open- 
circuit voltage Vo , and maximum enabling voltage V max • 
Hence, even if the cell parameter and open-circuit 
voltage Vo are instantaneously varied with respect to 

20 time (or continuously varied with respect to time), 
the output enabling power estimated value can be 
coincident with the real value. It is noted that an 
attention needs to be paid to develop the first-order 
lag of time constant T 3 in an apparent appearance on 

25 the real value of the open-circuit voltage and the 

terminal voltage (refer to right side second term of 
equation (18) of the cell model). 

[0059] In addition, in the input enabling power 

P ln of Fig. 91, a characteristic of reference (a dot- 
30 and-dash line) shown in Fig. 91 indicates a value 

calculated using the open-circuit voltage estimated 
value. As shown in Figs. 9A through 9 J, the input 
enabling power estimated value (dot-and-dash line) 
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calculated using the open-circuit voltage estimated 
value is larger than the real value of the input 
enabling power. This is caused by the fact that the 
apparent open-circuit voltage is larger than the real 
5 value of the open-circuit voltage Vo and is nearer to 
maximum enabling voltage V raax . In details, in a case 
where the input (charge) is carried out to the cell 
using the dot-and-dash lined input enabling power 
estimated value, there is a possibility that, in this 

10 case, the input enabling power real value breaks 

through maximum enabling voltage V max of the cell and 
the cell is deteriorated due to the overcharge. 
However, in the second embodiment, from the estimated 
coefficient parameters and open-circuit voltage Vo , 

15 the apparent open-circuit voltage Vo/C(s) 

(corresponds to AVo' in equation (35)) is calculated. 
Using one of Vo and Vo/C(s) which is nearer to 
maximum enabling voltage V max , the estimated 
parameters, and maximum enabling voltage V max , input 

20 enabling power P in is estimated. Hence, in the case 
of Figs. 9A through 9J, the input enabling power 
estimated value (solid line) is calculated using 
apparent open-circuit voltage Vo/C(s) nearer to 
maximum enabling power estimated value V max (solid 

25 line). Thus, the input enabling power estimated 

value is sufficiently coincident with the real value 
thereof and there is no possibility that the input 
enabling power estimated value breaks through the 
maximum enabling power voltage of the cell. 

30 [0060] On the other hand, in the column of output 

enabling power Pout in Fig. 9J, the characteristic 
described as the reference (dot-and-dash line) 
indicates a value calculated using the apparent open- 
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circuit voltage estimated value. As shown in Figs. 9J, 
the output enabling power estimated value (a dot-and- 
dash line) calculated using the apparent open-circuit 
voltage estimated value is larger than the real value 
5 of the output enabling power. This is caused by the 
fact that the open-circuit voltage estimated value is 
smaller than the apparent open-circuit voltage and is 
nearer to minimum enabling voltage V min . In details, 
in a case where the cell is outputted (discharged) 

10 using the dot-and-dash lined output enabling power 

estimated value, this output enabling power estimated 
value breaks through minimum enabling voltage V m in so 
that there is a possibility that the cell is 
deteriorated due to an over-discharge. However, in 

15 the second embodiment, the apparent open-circuit 
voltage Vo/C(s) is calculated from the estimated 
coefficient parameters and open circuit enabling 
voltage V rain . Using one of Vo and Vo/C(s) which is 
nearer to minimum enabling voltage V min , the estimated 

20 coefficient parameters, and minimum enabling voltage 
V m ±n. output enabling power P ou t is estimated. Hence, 
in the case of Figs. 9A through 9J, the output 
enabling power estimated value (solid line) is 
calculated using the open-circuit voltage Vo which is 

25 nearer to minimum enabling power V m±n . Hence, the 
estimated value of the output enabling power is 
sufficiently coincident with the real value and there 
is no possibility that the output enabling power 
estimated value breaks through minimum enabling 

30 voltage V min of the cell. It is noted that the relay 

described at steps S20A and S20B shown in Figs. 6 and 
7 corresponds to relay 70 shown in Fig. 2. 
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[0061] The entire contents of a Japanese Patent 

Application No. 2003-054035 (filed in Japan on 
February 28, 2003) are herein incorporated by 
reference. The scope of the invention is defined with 
reference to the following claims. 



